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Abstract  of  Thesis: 

Nanowire  metamaterials  are  a  class  of  composite  photonic  media  formed  by  an  array  of 
aligned  plasmonic  nano  wires  embedded  in  a  dielectric  matrix.  Depending  on  exact  composition, 
geometry,  and  excitation  wavelength,  nanowire  structures  are  known  to  exhibit  elliptical, 
hyperbolic,  or  epsilon-near-zero  (ENZ)  responses.  In  the  ENZ  regime,  optical  response  of  the 
composite  becomes  strongly  nonlocal.  Excitation  of  an  additional  wave,  caused  by  nonlocality, 
has  been  experimentally  demonstrated  in  nanowire-based  metamaterials.  In  this  thesis,  a 
computational  study  of  the  nonlocal  optical  response  in  plasmonic  nanowire  arrays  has  been 
conducted  to  better  understand  such  materials.  The  results  of  this  computational  study  were  used 
to  develop  an  analytical  technique  that  provides  an  adequate  description  of  the  optical  response 
of  wire  based  metamaterials.  This  fonnalism  combines  the  local  and  nonlocal  effective-medium 
theories  often  used  to  describe  the  optics  of  nanowire  composites.  It  provides  insight  into  the 
origin  of  the  additional  wave  and  allows  implementation  of  additional  boundary  conditions. 
This  approach  can  be  straightforwardly  extended  to  describe  the  optics  for  numerious  plasmonic 


structures. 
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1  Introduction 

Over  the  past  several  years,  there  has  been  an  increased  interest  from  the  scientific  community  to 
further  study  and  better  understand  metamaterials.  Metamaterials  are  a  class  of  composite 
materials  artificially  constructed  to  exhibit  exceptional  properties  that  are  not  readily  found  in 
nature.  There  are  numerous  applications  in  modem  optics  which  can  be  realized  through  the 
study  and  understanding  of  light  propagation  within  metamaterials.  These  range  from  imaging, 
sensing,  and  security,  to  solar  power,  optical  information  processing,  and  photonic  circuits  [1-9] 
One  such  metamaterial  that  has  recently  attracted  significant  attention  include  nanowire- 
based  composites,  a  class  of  materials  formed  by  an  array  of  aligned  plasmonic  nanowires 
embedded  in  a  dielectric  substrate  (Fig  1.1).  These  wire  materials  are  a  sub-class  of  uniaxial 
metamaterials  that  have  a  homogeneous  internal  structure  along  one  pre-selected  direction.  Due 
to  relatively  low  loss  and  ease  of  fabrication,  nanowire  composites  have  found  a  multitude  of 
applications  in  subwavelength  imaging,  biosensing,  acousto-optics,  and  ultrafast  all-optical 
processing,  spanning  visible  to  TFIz  frequencies  [10-15]. 


«z 


x 


Figure  1.1:  Schematic  geometry  of  wire  array. 
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Nanowire  structures  can  be  successfully  fabricated  using  electrochemistry  [7].  This  process 
generally  results  in  nanowire  composites  consisting  of  anodized  alumina  substrates  with  pre¬ 
etched  holes  filled  with  noble  metals,  such  as  gold  (Au)  [16].  For  these  types  of  materials,  the 
radius  r  for  the  individual  nanowires  can  be  controlled  between  10  and  100  nrn  with  a  separation 
distance  a  between  50  and  ~300«/?z. 

Not  only  do  nanowire-based  composites  fonn  the  platform  for  practical  applications, 
these  materials  have  also  gained  attention  due  to  their  unusual  and  counter-intuitive  optical 
properties.  These  behaviors  include  negative  refraction,  subwavelength  confinement  of  optical 
radiation,  and  modulation  of  photonic  density  of  states  [17-21].  The  unusual  optical  phenomena 
inside  these  nanowire-based  structures  can  often  be  accredited  to  extremely  large  or  vanishingly 
small  values  of  the  effective  refractive  index.  It  can  be  shown  that  at  these  extreme  values,  the 
validity  of  conventional  effective  medium  theories  (EMTs)  tends  to  be  inadequate  and  the 
effective  pennittivity  of  the  metamaterial  becomes  dependent  not  only  on  the  frequency  of  light, 
but  also  on  its  wavevector  [22-24].  This  becomes  obvious  in  the  epsilon-near-zero  (ENZ) 
regime,  where  the  contribution  of  spatial  dispersion  becomes  increasingly  important.  The 
presence  of  these  strong  non-local  terms  drastically  changes  the  optical  properties  of  the  material 
which  leads  to  the  appearance  of  an  additional  transverse  or  longitudinal  wave  [16]. 

The  broad  range  of  novel  fundamental  optical  phenomena  and  prospective  applications 
motivate  this  study  of  optical  properties  of  wire  materials.  The  rest  of  the  thesis  is  organized  as 
follows.  First,  a  description  of  effective  medium  theory  is  presented  along  with  a  discussion  of 
its  applicability  in  describing  wire  based  materials.  Second,  a  computational  and  numerical 
study  of  the  nonlocal  optical  response  in  plasmonic  nanowire  arrays  is  conducted  and  analyzed. 
This  acts  as  the  basis  of  comparison  for  all  other  numerical  and  analytical  calculations  to  follow. 
Third,  a  complete  investigation  of  the  numerical  results  and  a  comparison  with  the  current 
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models  is  conducted.  Finally,  an  analytical  technique  that  provides  an  adequate  description  of 
electromagnetism  in  wire-based  metamaterials  that  takes  into  account  nonlocal  optical  response 
originating  from  the  homogenization  procedure  is  obtained.  This  approach  can  be 
straightforwardly  extended  to  describe  optics  of  coaxial-cable-like  media  [25-27]  and  numerous 
other  uniaxial  composites. 
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2  Effective  Medium  Theory 

We  begin  with  one  of  the  most  common  approaches  in  analyzing  complex  optical  systems  with 
subwavelength  components.  This  is  known  as  effective  medium  theory.  In  this  technique,  the 
properties  of  composite  material  are  related  to  the  optical  properties  of  its  constituents.  The  small 
(in  comparison  to  wavelength)  size  of  the  components  of  the  material  makes  it  possible  to  solve 
Maxwell’s  equations  in  the  quasistatic  limit,  where  electric  and  magnetic  responses  of  the  system 
are  decoupled  from  each  other.  The  permittivity  of  (uniaxial)  wire  material  and  dispersion  of  the 
waves  propagating  in  this  medium  are  derived  below. 

2.1  Constitutive  relationships  in  the  quasistatic  limit 

We  assume  that  the  nanowire  system  operates  in  the  effective-medium  regime,  which  means  that 
its  inter-wire  separation  a  «  A0  with  A0  being  the  free  space  wavelength.  In  addition,  we  follow 
the  original  work  [28]  and  assume  that  surface  concentration  of  wires  is  relatively  small, 
p  —  nR2/a2  «  1,  where  R  is  the  radius  of  the  wire.  The  optical  response  of  the  nanowire 
materials  resembles  a  uniaxial  media  with  its  optical  axis  in  the  direction  of  the  wires.  As  such, 
the  permittivity  of  metamaterial  can  be  characterized  by  a  diagonal  tensor  with  components 

^ x,y  0 

0  £ x,y 

0  0 

Using  the  Generalized  Maxwell-Garnett  formalism  the  boundary  conditions  of  Maxwell’s 
equations  require  that  the  Ez  component  be  constant  across  the  unit  cell  Fig.  2.1a.  An  individual 
wire  excited  by  the  homogeneous  electric  field  in  the  z-direction  may  be  written  as 


: 

j  p-i-i] 
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Figure  2.1:  Dielectric  wire  in  a  constant  external  electric  field  (a)  electric  field  in  z-direction  and  (b)  in-plane 
electric  field  in  x-direction. 


pmg  _  mg 
cz  ez 


[2.1.2] 


The  in-plane  components  of  the  electric  field  (Fig.  2.1b)  can  be  derived  beginning  with  the 
solution  of  the  Laplace  equation  for  the  potential.  The  general  form  of  the  potential  for  a 
homogeneous  material  in  cylindrical  coordinates  is 


A  [2.1.3] 

d>  =  (A0  +  B o  lnr)  +  ^  (A- mrm  +  Bmr  m)(Cm  cos (m0)  +  Dm  sin (m0)) 

m=  1 


Where  the  coefficients  of  linear  expansion  AL,  Bit  CL  and  DL  are  determined  by  the  requirement 
that  the  potential  is  finite  at  r  =  0,  that  the  field  away  from  the  wire  is  converging  to  the  field  of 
the  plane  wave  »  R)  —  — e™5r  cos  c/>J,  and  by  the  boundary  conditions  requiring  the 
continuity  of  Et(j  —  R)  and  D±(r  —  R).  Therefore,  inside  the  wire  r  <  R 

^  [2.1.4] 

d>j  =  A  +  y  rm(Cm  cos (m<p)  +  Dm  sin (m<p))  —  A  +  rC1  cos  <p 

m= 1 
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d<$>i  1 

Et  —  — V<t>j  =  — - — r - -—<6  =  C1(—  cos  <p  r  +  sin  (p  (f>) 

or  r  dtp 


while  outside  the  wire  r  >  R 


[2.1.5] 


<A>h  —  A  +  ( A yr  +  Byr  x)  cos  0 


[2.1.6] 


[2.1.7] 


It  is  instructive  to  rewrite  the  expressions  above  in  Cartesian  coordinates, 


Ei  = 


2eh 


£i  +  £h 


[2.1.8] 


^1  +R2 


f n  -  gf  y2  -  *2  \ 

6i  +eh  ( x 2  +  y2Y) 


x  — 


2R2xy  eh  -  et  _ 
(x2  +y2)2£i  +£hy_ 


[2.1.9] 


Note  that  in  the  considered  case  of  the  excitation  field  is  directed  along  the  x  axis,  and  the  y 
component  of  the  excited  field  is  anti-symmetric  with  respect  to  x,  y.  Therefore,  the  only  relevant 
component  that  does  not  vanish  when  averaged  over  the  unit  cell  of  the  field  becomes  (see 
however,  [30]) 


E?a  =  e7  x  { 


£i  +  £n 


,r  <  R 


1  +  R‘ 


eh 


y2  —  x2 
2)2 


et  +  eh  (x2  +  y2) 


,r 


[2.1.10] 


Here  £[  and  eh  being  the  permitivitties  of  the  wire  inclusion  and  of  the  host  material  respectively, 
and  the  parameters  e™5and  e™5  are  the  field  amplitudes.  In  the  limit  of  small  wire  concentration, 
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interaction  between  the  wires  can  be  neglected.  Taking  the  average  of  the  fields  over  the  unit  cell 
for  the  j- th  component  yields  the  effective  permittivity  ejl'g  —  ( e (x,  y) E™a (x,y))/{ E™9  (x,y)): 


mg  _ 

tZ 


Ve^a  +  (1  -  p)eht?a 

ea 


=  P^i  +  (1  -p)eft 


[2.1.11] 


mg 

Lx,y 


P^i 


2 £h  mg 

ei  +  eh  x 


2  Eh 
ei  + 


mg 


+  (1  -pKe™5 

+  d  ~pK3 


2pejeh  +  (1  ~  pK(fi  +  fh) 

2peh  +  (1  -p)(ef  +  eh)  [2.1.12] 


The  derivations  outlined  above  was  originally  considered  in  the  context  of  the  permittivity  for 
nanosphere-based  composites[28],  often  referred  to  as  Maxwell  Garnett  formalism. 

2.2  Dispersion  Relationships 

In  this  section  let  us  derive  the  dispersion  relationships  governing  light  propagation  in 
homogeneous  uniaxial  materials,  and  therefore  governing  the  propagation  of  light  in  nanowire 
metamaterials  in  the  limit  of  effective  medium  theory  (EMT).  We  begin  with  Maxwell’s 
equations 
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V  •  D  —  Anp 


V  •  B  =  0 


V 


x  E  =  - 


ldB 
c  dt 


If  dD N 

V  x  H  =  -  [AnJ  +  — 
c\  dt  , 


For  an  anisotropic  composite  the  dielectric  tensor  must  be  diagonal  and  was 
Eqn.(2. 1.1).  The  propagating  plane  waves  can  be  written  as 


E  =  £’0exp(ifc  •  r  —  ia)t ) 


H  —  H0exp(ik  ■  f  —  icdt) 


By  implementing  Eqn.(2.2. 1)  with  Eqn.(2.2.5)  we  get 


k  ■  eE0  —  0 


[2.2.1] 

[2.2.2] 

[2.2.3] 

[2.2.4] 

shown  in 

[2.2.5] 

[2.2.6] 

[2.2.7] 


and  Eqn.(2.2.3)  with  Eqns.(2.2.5)  and  (2.2.6)  we  can  write  the  magnetic  field  in  anisotropic  non¬ 
magnetic  medium  as 
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C 


[2.2.8] 


<D 


Finally  Eqn.(2.2.4)  on  Eqn.(2.2.6)  and  the  relation  e0/i0c2  =  1  we  obtain 


Eqns.(2.2.7)  and  (2.2.9)  make  up  Maxwell’s  equations  for  a  nonmagnetic  anisotropic  material. 

Note  that  D  1  k,  D  1  H,  H  1  E,  H  1  k.  It  was  shown  above  that  the  vectors  D  and  E 
are  not  necessarily  co-aligned  in  anisotropic  media,  so  E  and  k  are  not  necessarily  orthogonal  to 
each  other,  see  Fig.  2.2c. 


B,H 


Figure  2.2:  (b)  TE-polarized  ordinary  waves  are  identical  to  isotropic  systems,  (c)  TM-polarized  extraordinary 
waves  exhibit  hyperbolic  dispersion  and  have  no  diffraction  limit  when  exy  ■  ez  <  0  . 

One  can  obtain  from  Eqns.(2.2.8)  and  (2.2.9)  the  following  relationship  where  any  plane  wave 
that  satisfies  Maxwell’s  equation  must  also  satisfy  the  expression 


k2E o  —  •  E0 ) 


[2.2.10] 
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with  k2  —  k  ■  k.  In  component  form,  where  indices  i,j  correspond  to  Cartesian  components, 
Eqn.(2. 1.10)  can  be  written  as 


nr 


Sijk2  -  kikj  -  6tj  —  1  E0i  =  0 


[2.2.11] 


In  order  for  the  electric  field  to  have  a  non-zero  amplitude  the  dispersion  relation  for  a  wave 
propagating  in  an  anisotropic  material  reduces  to  the  well-known  expression 


(0‘ 


Sijk  kikj  6ij  ^ 


=  0 


[2.2.12] 


with  k  and  co  being  the  wave  vector  and  angular  frequency  of  the  plane  wave,  respectively,  c 
being  the  speed  of  light  in  vacuum,  e  is  generally  the  nonlocal  dielectric  pennittivity  tensor  of 
the  metamaterial,  and  is  the  Kroenecker  delta  symbol.  When  the  pennittivity  does  not 
depend  on  wavevector  (also  known  as  local  regime)  e  —  e™3 .  Expanding  Eq.(2.2. 12) 


k2  +  k2x-e™^ 


~kxky 


kXky 


nr 


-kxkz 


k  z  4-  k  2  -  F  y  _ 

a z  cx,y  ^2 


kxkz 


l^y 


k2  +  k2  -  F"ly  — 

n,x  i  a .y  Cz 


-  0 


[2.2.13] 


and  solving  for  the  dispersion  allows  two  families  of  waves.  One  of  these  is  known  as  ordinary, 
or  transverse-electric  wave,  which  is  not  affected  by  material  anisotropy.  While  the  second  is 
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often  called  the  extraordinary,  or  transverse-magnetic  wave  (see  Fig.  2.2).  The  dispersion  of 
these  waves  is  characterized  as: 


kx  +  ky  +  k2 


[2.2.14] 


k2  k2  +  ky  co2 

mg  d  mg  ~  [2.2.15] 

tx,y  tz 


By  adjusting  composition  of  the  metamaterial  and  operating  wavelength,  the  optical  response  of 
the  composite  can  be  controlled  between  all-dielectric  elliptic  (e™9  >  0,  e™9  >  0)  epsilon-near- 
zero  (ENZ,  e ™9  ~  0)  and  hyperbolic,  which  consists  of  (e™9  >  0,  e ™9  <  0)  type  I  and  ( e ™9  < 
0,  e ™9  >  0)  type  II  regimes.  In  the  two  latter  regimes,  metamaterial  supports  optical  waves  with 
either  small  or  large  effective  modal  index  and  motivate  a  number  of  potential  applications  in 
molding  of  light  [31-33],  cloaking  [4,5],  and  subwavelength  light  manipulation^  1-13,17]. 


Figure  2.3:  £-space  topology.  The  dispersion  for  (a)  and  isoptropic  dielectric  is  a  sphere,  and  for  (b)  extraordinary 
waves  in  a  uniaxial  medium  with  extreme  anisotropy  it  is  a  hyperboloid  (type  I).  (c)  Hyperboloid  of  a  type  II 
metamaterial. 


To  reiterate,  a  unique  feature  of  hyperbolic  metamaterials  is  evident  from  Fig.  2.3,  which 
shows  that  the  wavevector  can  exceed  the  free  space  wave  vector  (/c0  =  m/c).  Hyperbolic 
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metamaterials  support  a  large  number  of  electromagnetic  states.  An  intuitive  counting  procedure 
in  A-spacc  consists  of  calculating  the  volume  between  the  dispersion  contours  at  <u(/c)  and 
6j(/c)  +  Am.  In  a  vacuum  it  is  the  equivalent  to  an  infinitesimally  thin  spherical  shell  in  A-space. 
For  a  hyperbolic  metamaterial  the  dispersion  modes  have  been  shown  to  lie  within  the  surface  of 
a  hyperboidal  shell  which  has  an  infinite  volume.  In  the  effective  medium  limit  this  leads  to  a 
broadband  singularity  in  the  density  of  states  [34],  There  is  no  upper  cut-off  to  the  wave  number 
in  an  ideal  hyperbolic  medium  and  an  infinite  number  of  A  modes  are  supported.  This 
phenomenon  is  the  drive  for  further  investigation  and  has  vast  potential  in  application  for 
hyperbolic  metamaterials. 

2.3  Nonlocality 

For  most  materials  the  dielectric  permittivity  tensor  e  can  be  treated  as  a  function  of  frequency 
alone.  However,  with  the  material  sizes  becoming  smaller  and  smaller  in  length  scales  the  finite 
wavelength  of  electrons  can  no  longer  be  ignored.  Following  Landau  and  Lifshits  [51]  (see  also 
[22,  35]),  it  is  more  correct  and  appropriate  to  use  an  approach  based  on  taking  spatial  dispersion 
into  account.  An  electric  field  in  one  place  can  produce  a  polarization  in  the  near  vicinity, 

P(r)  =  jx(r~r')E(r')dr'  I23'1! 


By  converting  Eq.(2.3.1)  in  reciprocal  space  we  obtain  the  following  fonn 

P(fc)  =  x(k)E(k)  [2-3.2] 


where  polarizability  (and,  consequently,  pennittivity)  depends  on  the  wavevector  k.  Note  that  if 
the  material  response  is  local  [x(r )  oc  8(r  —  r')],  pennittivity  in  wavevector-domain  does  not 
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explicitly  depend  on  k.  In  a  majority  of  realistic  materials  nonlocality,  if  present,  is  weak,  and 
thus  it  can  be  considered  to  be  a  correction  to  frequency-dependent  response.  If  we  take  a  Taylor 
expansion  of  e(k),  in  centrosymmetic  materials,  the  lowest  non-vanishing  term  is  proportional  to 
k2. 


3  3 

eij(k)  =  ej(<u)  ^  cCijrniWkikm 

1  =  1  771=1 


[2.3.3] 


where  is  a  4th  rank  tensor. 

The  spatial  dispersion  appears  in  addition  to  the  more  familiar  temporal,  or  frequency, 
dispersion,  which  is  reflected  in  the  dependence  of  the  tensor  on  the  frequency,  a).  The  spatial 
dispersion  effects  are  usually  much  weaker  than  those  arising  from  the  frequency  dispersion,  but 
they  can  lead  to  qualitatively  new  phenomena,  by  substituting  Eq.  (2.3.3)  into  Eq.  (2.2.12)  an 
additional  electromagnetic  wave  arises.  These  phenomena  become  increasingly  important  in  the 
regime  where  permittivity  of  material  approaches  zero  (often  known  as  epsilon-near-zero,  ENZ 
regime)  since  in  this  limit  spatial  dispersion,  however  small,  dominates  the  material  response. 

Similar  to  homogeneous  materials,  it  was  shown  that  for  nanowire  meta-materials 
nonlocality  needs  to  be  taken  into  consideration  at  least  in  the  ENZ  regime[16].  Here  the  spatial 
dispersion  changes  the  optical  properties  of  the  structure  which  leads  to  the  excitation  of  an 
additional  transverse-magnetic  wave  that  is  not  described  by  local  EMTs. 
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3  Numerical  Calculations 

For  simplicity,  the  frequency  of  electromagnetic  excitations  and  the  unit  cell  parameters 
of  the  system  are  fixed  and  only  the  permittivity  of  the  wire  inclusions  or  substrate  host  are 
varied.  The  parameters  primarily  used  throughout  this  work  are  R  —  20  nm,  a  —  100  nm, 
eh  —  1,  L  —  1  ixm,  A0  —  1.5 [im  (see  Fig.  3.1),  which  are  typical  for  composites  fabricated  with 
anodized  alumina  templates  [7,36].  Permittivity  of  the  wire  was  varied  between  -10  and  -3.  This 
ensures  that  the  parameters  of  the  composite  cover  the  important  transition  between  elliptical 
response  (fj  ^  —7),  epsilon-near-zero  (ENZ)  regime  ~  —7),  and  hyperbolic  response 
(<fi  <  -7)  [see  Fig. 3. 2]. 

Two  separate  numerical  techniques  for  solving  Maxwell’s  equations  in  periodic 
geometries  were  used,  finite  elements  technique  (FEM)  and  rigorous  coupled  wave  analysis 
(RCWA). 


.y 

i  i 


Figure  3.1:  Schematic  geometry  and  a  unit  cell  of  a  nanowire  composite. 

First,  a  FEM  solver  [37]  was  used  to  perfonn  numerical  analysis  of  the  eigenmodes  of  the 
periodic  nanowire  composites.  A  single  unit  cell  was  constructed  as  shown  in  (Fig.  3.1).  To 
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mimic  the  response  of  the  infinite  periodic  wire  medium,  we  used  Floquet  periodicity  for  the 


boundary  conditions. 


Figure  3.2:  Permittivity  derived  from  local  EMT,  e (dashed  line)  and  e™a  (solid  line).  The  ENZ  regime  for  the 
particular  geometry  under  investigation  is  at  f;  =  —  7  and  illustrated  by  the  cross  in  (a). 


As  with  any  periodic  composite  (photonic  crystal),  nanowire  arrays  support  a  large 
number  of  modes  that  essentially  resemble  the  series  of  Bloch  waves.  It  has  been  previously 
shown  [16]  that  in  the  limit  of  small  lattice  size,  depending  on  the  material  parameters,  optical 
performance  of  nano  wire  composites  is  dominated  by  either  two  or  three  modes.  The  properties 
of  one  of  these  waves  resemble  the  properties  of  TE-polarized  mode  in  homogeneous 
composites,  while  the  unit-cell-averaged  polarization  of  the  two  other  waves  suggest  their 
transverse-magnetic  (TM)-polarized  character,  Fig.  3.3. 

Due  to  the  existence  of  multiple  beams  with  identical  polarization  it  is  reasonable  to 
assume  that  the  behavior  of  the  two  average-TM-polarized  waves  in  wire  composites  can  be 
described  by  nonlocal  (wavevector-dependent)  effective  medium  theory.  The  goal  of  this  thesis 
is  to  develop  such  a  description. 
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Figure  3.3:  (a-c)  Electric  field  and  (d-f)  Magnetic  field  profiles  produced  using  FEM  calculations.  (a,d)  Represents 


the  TE-mode,  (b,e)  the  TM  or  extraordinary  -mode  and  (c,f)  the  additional  TM-mode  which  has  been  predicted  and 


observed  in  previous  experimental  results. 


In  this  model,  the  two  TM  waves  identified  in  the  numerical  simulations  represent  the 


main  and  additional  TM  waves  that  are  supported  by  the  nanowire  metamaterial.  The  main  wave 


has  a  relatively  smooth  Ez  profile  (Fig.  3.3  b,e)  in  comparison  to  the  additional  wave  which  has 


strong  non-uniform  field  profile  (Fig.  3.3  c,f). 


3.1  Two  Dimensional  Rigorous  Coupled  Wave  Analysis  (2DRCWA) 

To  provide  independent  verification  of  FEM-based  solutions  of  Maxwell  equations,  the  2D 
analog  of  RCWA  [38]  is  developed  and  implemented.  In  this  formalism,  that  takes  advantage  of 
2D  periodicity  of  material  permittivity,  the  latter  is  represented  using  2D  discrete  Fourier 
transform: 
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e 


(x,y)  =  V 


- i(mqxx+nqyy ) 


m,n 


[3-1.1] 


where  qa  —  —  plays  the  role  of  the  Bloch-vector,  and  subscript  a  defines  the  in-plane  coordinate 


(. x  or  y).  Using  the  Bloch  theorem,  we  represent  the  fields  in  the  unit  cell  by 


E,H  -  '^\e,h)rnneikrn’n'r~i(0t 

m,n 


with 


kjn,n  —  (jix  -|-  rnq^x  +  (/cy  -I-  nqy)y  +  kzz 


[3.1.2] 


[3.1.3] 


Explicit  substitutions  of  Eqs.(3.1.1  and  3.1.3)  into  Eq.(2.2.12),  yields  the  dispersion  relation 


a,/3 


i  r  i2 

’  |  $a-m,p-n  ^2 ~^a-m,[3-n 


=  0 


[3.3.4] 


Replacing  ( m,n )  — >  and  (a,/?)  ->  rj  transforms  the  mxn  and  ax/?  matrices  into  a  diagonal 
matrix,  where  the  indices  and  17  are  one-dimensional  indices  that  now  span  two-dimensional 
space.  With  this  transformation  and  the  use  of  Maxwell’s  equations  we  obtain 


0  X  hfr  —  — 


CO 


[3.3.5] 


/r  x  ec  —  —  hr 

s  c  s 


[3.3.6] 
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The  Eqns.(3.3.5)  and  (3.3.6)  now  fonn  a  basis  where  the  eigenvalue  kz  and  eigenvectors  e,h  for 
each  mode  can  be  solved.  The  FEM  simulations  and  2D  RCWA  yield  comparable  results  (Fig. 
3.4). 


Figure  3.4  Dispersion  of  the  modes  guided  by  the  nanowire  metamaterial,  calculated  with  full-wave 
solutions  of  Maxwell  equations  (lines)  and  2DRCWA  (symbols) 

3.2  Comparison  of  Numerical  Results  with  EMT 

To  first  explore  the  validity  and  applicability  of  EMTs  the  dispersion  of  the  TE -polarized 
and  TM-polarized  waves  are  extracted  from  numerical  calculations  and  compared  with  the 
dispersion  predicted  by  Eqs.(2.2.14)  and  (2.2.15).  One  way  that  the  parameter  exy  and  ez  can  be 
determined  numerically  is  by  minimizing  the  mean-square  deviation, 
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A= 'Ei\k-zi  —  kz\2  min.  [3.2.1] 

where  kz  is  calculated  using  Eq.(2.2.14)  or  (2.2.15),  while  kz[  is  taken  from  the  numerical 
solutions  to  Maxwell’s  equations  (Fig.  3.4).  Eqn.(3.2.1)  can  be  solved  analytically  from  the  TE- 
polarized  wave  yielding 


i=  1 


[3.2.2] 


Similarly,  the  parameter  ez  can  be  detennined  by  minimizing  the  deviation  of  calculated 
dispersion  of  main  TM-polarized  wave  from  its  expected  elliptical  (or  hyperbolic)  behavior.  This 
procedure  yields 
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-LS 


-1 


i= i  e 


LS 

x,y 


0) 


™ _ L-2 


[3.2.3] 


Once  the  effective  permittivity  is  fit  to  the  numerical  solutions  of  Maxwell’s  equations,  a 
comparison  can  be  made  to  the  analytical  predictions  of  EMT.  This  comparison  between  the 
mean-square  deviation  and  analytical  effective  medium  parameters  is  shown  in  Fig.3.5.  It  is  seen 
that  exy  is  adequately  described  by  Eq.  (2.1.12).  However,  Eq.  (2.1.11)  deviates  from  the  fit 
value  and  this  deviation  is  particularly  strong  at  the  ENZ  ( ez  — >  0)  regime.  In  particular, 
numerically-derived  ez  seems  to  depend  on  the  kx  component  of  pennittivity;  with  Eq.(2.1.11) 
less  accurately  describing  optics  of  nanowire  systems  for  larger  values  of  kx.  This  is  validated 
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by  reducing  the  limit  of  the  maximum  kx  term  used  in  the  summation  to  determine  the  numerical 
permittivities,  Eqs.(3.2.2)  and  (3.2.3).  As  this  tenn  is  reduced  it  can  be  seen  from  Fig.  3.5  that 
the  permittivity  ez  (triangles)  is  now  in  better  agreement  with  EMT  (solid  line). 

(a)  (b) 


Figure  3.5  Permittivity  derived  from  local  EMT  (lines)  compared  with  the  permittivity  calculated  using  least  square 
approximation  (symbols)  eqns.  3.2.2  and  3.2.3.  The  transverse  permittivity  is  in  good  agreement  for  all  values  of  kx 
however  the  permittivity  parallel  to  the  wires  is  not  especially  near  ENZ  (squares).  The  agreement  improves  for 
smaller  maximum  value  of  kx  (triangles). 


When  ez  is  sufficiently  different  from  zero,  however,  there  seems  to  be  reasonable  agreement 
between  the  predictions  of  local  EMT  [6,30,39,40]  and  properties  of  the  main  TM-polarized 
wave  in  wire  composites,  in  particular  in  the  limit  of  relatively  small  kx  (Figs.  3. 5-3. 7). 

To  further  investigate  the  relationship  between  the  numerical  and  the  local  dispersion, 
Eqs.(2.1.11)  and  (2.1.12),  are  plotted.  It  can  be  seen  from  Fig.  3.6  that  the  local  theories 
inadequately  describe  the  dispersion  relationship  for  the  transverse-magnetic  polarized  wave. 
Finally  the  numerical  permittivity  obtained  by  averaging  the  wave  profiles  from  the  FEM  results 
and  the  analytical  EMT  permittivity  as  a  function  of  kx  are  compared.  Here  (Fig.  3.7)  the 
deviation  at  large  values  of  kx  become  obvious  for  the  TM-mode. 
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Figure  3.6  A  comparison  between  local  EMT  dispersion  (dashed  lines)  and  FEM  calculations,  (symbols).  (a,b)  TE- 
mode  comparison  which  is  in  good  agreements  (note  the  axis  is  shifted  by  0.1  to  show  distinction  between  each 
dispersion  curve).  (c,d)  TM-mode  comparison,  here  it  can  be  clearly  shown  that  there  is  a  deviation  for  large  values 
of  kr. 
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Figure  3.7  A  comparison  between  local  EMT  permittivity  (lines)  and  FEM  calculations,  (symbols),  (a)  TE-mode 
comparison  which  is  in  good  agreements,  (b)  TM-mode  comparison,  here  it  can  be  clearly  shown  that  there  is  a 
deviation  for  large  values  of  kx . 
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To  summarize,  local  effective  medium  theory  falls  short  in  several  major  areas.  First,  it 
fails  to  accurately  describe  the  fundamental  optical  properties  for  uniaxial  materials  at 
relatively  large  (in  comparison  with  oo/c )  values  of  kx,  and  near  the  ENZ  regime.  Second, 
and  more  importantly,  it  does  not  predict  the  additional  TM-mode.  In  the  next  sections  we 
summarize  previous  attempts  to  develop  nonlocal  effective  medium  theories  that  attempt  to 
describe  optical  properties  of  wire  composites. 

3.3  Comparison  with  earlier  nonlocal  effective  medium  theories 

The  problem  of  light  interaction  with  nonlocal  wire  media  has  been  considered  before.  A 
realistic  wire  composite  in  the  effective  medium  regime  supports  three  electromagnetic  modes. 
With  few  exceptions,  the  majority  of  the  studies  that  focus  on  high-frequency  response  of  wire 
composites  [6,30,39-41]  predict  a  single  extraordinary  and  a  single  ordinary  wave  at  every 
frequency.  Most  of  the  studies  that  do  predict  the  existence  of  additional  electromagnetic  waves 
[24,42,43]  have  focused  on  optics  for  highly  conductive  (PEC-like)  wires.  Refs.  [24,42,43] 
presents  an  attempt  to  generalize  the  developed  formalism  to  the  case  of  plasmonic  media.  Figs. 
3.8  and  3.9  present  a  comparison  of  the  formalism  from  [24,43]  and  full  wave  numerical 
solutions  of  Maxwell’s  equations. 

Fig.  3.8  clearly  demonstrates  that  the  approach  developed  in  [24,43]  severely 
underestimates  the  effective  modal  index  of  the  waves  propagating  in  plasmonic  wire  media 
(similar  phenomenon  can  also  be  seen  in  Fig.  3  of  [24]).  This  under-estimation  yields  significant 
errors  in  calculations  of  optical  properties  of  wire  composites,  seen  in  Fig.  5  of  [43]  that  can  be 
only  be  eliminated  by  using  a  heuristic  “correction  factor”. 
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kxc/w  kxc/w 


Figure  3.8:  Dispersion  of  the  modes  guided  by  the  nanowire  metamaterial,  calculated  with  lull-wave  solutions  of 
Maxwell  equations  (symbols)  with  earlier  approach,  designed  for  highly-conducting  wires  [24,43]  (dashed  lines) 

The  [24,43]  model  and  numerical  effective  pennittivity  for  increasing  kx  is  also  shown  to 
have  significant  errors.  As  kx  increases  the  variation  becomes  dramatic  and  the  need  for  a  more 
accurate  model  could  not  be  more  evident. 
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Figure  3.9:  Effective  permittivity  defined  as  ez  =  ( e(x,y)Ez(x,y))/(Ez(x,y )),  calculated  with  full-wave  solutions 
of  Maxwell  equations  (symbols)  with  earlier  approach,  designed  for  highly-conducting  wires  [24,43]  (dashed  lines). 
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4  Beyond  Effective  Medium  Theory:  Perturbation  approach 

One  major  problem  is  that  the  deviation  of  dispersion  of  the  main  TM-polarized  mode  in 
wire  composites  from  the  predictions  of  local  EMT  becomes  much  stronger  in  the  ENZ  regime. 
This  behavior  is  qualitatively  similar  to  behavior  predicted  (and  later  observed  at  ultra-low 
temperatures  [22,23,44])  in  nonlocal  (pennittivity  dependent  on  wavevector)  homogeneous 
materials.  Quantitatively,  (weak)  dependence  of  effective  pennittivity  on  wavevector  can  be 
considered  with  a  perturbation-type  conection  to  frequency-dependent  permittivity,  given  by 


kW 

£z  —  £ z  (<U)  T  2 


[4.1] 


with  \az\  «  1. 

The  nonlocality  of  the  z  component  of  the  pennittivity  only  affects  propagation  of  the 
TM-polarized  mode.  This  wavevector-dependent  term  changes  the  response  of  the  system  thus 
enabling  a  new  kind  of  plane  wave,  known  as  the  additional  wave.  By  substituting  Eq.  (4. 1)  into 
(2.2.13)  a  modified  dispersion  relationship  can  be  derived  that  takes  into  account  this  conection, 
and  is  given  by 


or 


O) 


[4.2] 


2azki  =  (ez  +  ^)  —  +  -  (ez  -  <f)2  —  +  4 Qc/ 


where  (  =  cczeXy 

In  the  limiting  case  where  az  — >  0,  one  of  the  roots  converges  to  Eq.  (2.2.15)  while  the 
other  root  goes  to  infinity.  Therefore,  the  additional  wave  is  always  present  in  nonlocal 
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metamaterials,  but  the  large  index  mismatch  may  prevent  its  coupling  with  free  plane  waves. 
The  ratio  between  the  parameter  az  and  the  minimal  value  of  ez(u>)  is  the  effectiveness  of 
coupling  incident  light  into  the  additional  wave.  Excitation  of  the  additional  wave  at  ENZ  in 
response  of  wire  composites  has  been  experimentally  demonstrated,  at  room  temperature,  in 
[16]. 

The  fitting  procedure  that  was  developed  in  Section  3.2  is  implemented  to  extract  the 
parameters  of  perturbation-based  nonlocal  dispersion  relation,  Eq.  (4.2).  Note  that  the  least 
square  fitting  procedure  becomes  nonlinear;  therefore,  empirical  solutions  for  the  four 
parameters,  5He(ez),  3m(ez),  5Re(az),  and  3m(az),  cannot  be  determined.  Instead  numerical 
values  of  these  four  parameters  were  simultaneously  determined  by  nonlinear  least  square  fit 
optimization[45-50] . 

After  a  solution  for  the  parameters  are  determined  the  nonlocal  dispersion  can  be 
calculated  and  compared  with  the  numerical  solutions  from  Maxwell’s  equations  Fig.  4.1.  It  can 
be  seen  in  the  figure  below  that  this  technique  appears  to  predict  both  waves  and  describes  the 
dispersion  with  excellent  accuracy. 
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Figure  4.1:  Dispersion  of  the  modes  guided  by  the  nanowire  metamaterial,  calculated  with  full-wave  solutions  of 
Maxwell  equations  (symbols)  with  nonlocal  dispersion  approach  (lines). 

However,  further  analysis  reveals  that  the  nonlocality  parameter  az,  as  determined  by 
nonlinear  optimization  procedure,  is  mode-dependent  (Fig.  4.2),  i.e.  it  has  different  values  for 
main  and  additional  TM-polarized  modes.  Thus,  it  is  obvious  that  perturbative  technique  Eq. 
(4.2)  does  not  provide  an  adequate  description  of  optical  properties  of  plasmonic  wire  arrays. 


27 


38 


(a)  (b) 


Figure  4.2:  (a)  Comparison  between  e™9  (line)  and  e™in  (symbols).  The  triangles  represent  the  TM  or  extraordinary 
mode  and  squares  the  additional  TM-mode.  (b)  The  nonlocal  fitting  parameter  (symbols)  for  varying  et,  here 
the  mode  dependency  can  be  clearly  seen. 


28 


39 


5  Analytical  Description  of  the  Nonlocal  Optical  Response 

In  this  section  an  analytical  description  for  the  optical  properties  of  nonlocal  nanowire  systems  is 
presented.  First,  light  propagation  in  the  direction  of  the  wire  is  considered  and  the  propagation 
constants  of  the  two  TM-polarized  waves  that  are  supported  by  the  composite  are  identified  in 
this  regime.  Lastly,  it  is  shown  how  to  calculate  the  dispersion  and  wave  profiles  for  off  angle 
light  propagation  and  compare  these  results  with  numerical  calculations. 

The  electromagnetic  field  within  the  unit  cell  (Fig.  3.1)  can  be  derived  from  first 
principles.  Within  the  nanowire  material  there  are  no  free  charges  and  Maxwell’s  equations  are 
used  to  obtain  the  generalized  wave  equation 


V2F 


=  o 

cz  dt2 


[5.1] 


where  F  =  E,  H.  Since  the  nano  wire  material  is  nonmagnetic,  without  loss  of  generality,  we  can 
set  the  permeability  to  unity.  The  wave  equation  then  reduces  to  the  following  fonn 


V2F 


—  —  =  0 

c2dt2 


[5.2] 


It  can  be  shown  (see  below)  that  for  a  wave  propagating  along  the  z  direction,  all  the 
transverse  components  of  F  can  be  expressed  in  terms  of  the  longitudinal  component  [51].  In 
Cylindrical  coordinates  the  differential  equation  for  this  remaining  component  takes  the  form 


V2Fz  -  k2Fz 


[5.3] 
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where  k2  —  eco2 /c2 .  The  solution  can  be  expressed  as  the  product  Fz  —  P(r)P(0)Z(z). 

Making  this  substitution,  the  differential  equation  can  then  be  written  as 

ii(r^4^  +  i^  +  fe2=  o.  ^ 

rR  dr  V  dry  r2P  dc/)2  Z  dz2 

From  the  above  equation  we  arrive  at  the  solution  Z(z)  =  z0elkzZ,  and  if  we  multiply  through  by 
r2  we  are  left  with 


r  d/dR\ 
R dr  V  dr  J 


1  dzP 


+  »^+<<k-kz)r2  =  0. 


[5.5] 


Notice  that  the  equation  is  now  completely  separable  by  variables.  Let’s  first  solve  for  P, 


1  d2P 
P  d02 


—  A2  —  const. 


[5.6] 


we  can  choose  A2  —  —m2  with  m  =  0,1,2,..  leading  to 


P(0)  =  Asin(rn<p)  +  Bcos(rn<p ). 


[5.7] 


The  differential  equation  can  now  be  written  in  terms  of  only  the  variable  r 


r  d 
R  dr 


—  m2  +  (k2  —  k2)r2 


n  d2R  dR  n  n  n 
—  r2  r  +  r  —  +  [(/c2  —  k2)r2  —  m2]R  =  0 
dr2  dr 


[5.8] 


Let’s  first  choose  k  =  yj ecj2 /c2  —  k2  so  the  equation  becomes 
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Here  boundary  conditions  must  be  applied  to  solve  for  the  amplitudes  of  the  individual 
cylindrical  waves. 

5.1  Calculation  of  Ez  and  Hz 

The  nanowire  metamaterial  can  be  separated  into  two  regions,  (i)  inside  and  (ii)  outside  the  wire 
material.  An  expression  for  the  fields  will  be  derived  for  both  by  applying  initial  boundary 
conditions  to  the  generalized  field  expression  previously  obtained  (Eq.  5.11). 

Inside  the  nanowire 

First  we  will  look  at  the  solution  inside  the  wire,  we  can  write  kl  —  eL(jj2 /c2  —  kz  and  as 
r  -»  0,  Em(/r,r)  00  thus  c2m  and  c4m  =  0.  Due  to  the  symmetric  nature  of  the  problem  it  is 
satisfactory  to  examine  the  fields  at  just  z  =  0, 
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Fl 

1  z 


V  JmOctnrXMn  sin(m0)  +  Bin  CO s(m0)) 

m 


[5.1-1] 


and  fz  oc  cos (m<p)  and  oc  sin(m0),  using  these  results  the  inner  wire  field  equations  can  be 
written  as 


Elz(r  <  fl)  =  ^  anJm(Kinr')  cos (m0) 

m 

Hlz(r  <  R)  —  ^  bmJmiKinr)  Sin (m0) 


[5.1.2] 


[5.1.3] 


Outside  the  nanowire 

Outside  the  wire  the  second  order  Bessel  function  does  not  vanish.  We  can  express  this  result  in 
tenns  of  Hankel  functions,  where  Kh  —  sj eha)2 /c2  —  kz 


Fn 

1  Z 


=  'Y(c?mH+(Khr)  +  c£mtfm(Khr))(2l^sin(m0)  +  cos (m0)) 


[5.1.4] 


where  H ^  —  Jm  +  iYm  and  Hm  —  Jm  —  iYm.  The  outer  wire  fields  can  be  written  as 


Elz(r  >  R)  =  ^\ainH^{Khr)  +  amHm(Khr))  cos (m0) 

m 

Hlz(r  >  R)  =  (Khr)  +  PmHm(.Khr))  sin(m0) 


[5.1.5] 


[5.1.6] 
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5.2  Calculation  of  Er ,  //r,  E and  H ^ 


To  obtain  the  fields  in  the  other  two  directions  we  must  begin  with  Maxwell’s  equations 

_  _  ldH  .  _  „  ed E  TT  .  .  .  .  .  _  .  ^(ldAz  dA <fc\  ~  fdAr 

V  x  E  = - —  and  V  x  H  —  Using  the  identity  V  x  A  =  r  ( - -)  +  p  ( - 

c  dt  c  dt  b  J  \r  dip  dz  J  ^  \  dz 

)  +  z  (“£:  (jAq')  —  ~~^j  we  can  write  down  the  following  relationships 


dAz 

dr 


dEr  c  | 

fldHz 

dH 0\ 

[5.2.1] 

dHr  . 

1  dEz 

dt  e ' 

\r  dp 

dz  J 

dt  _C' 

^  dz 

r  dp 

dEd 


c  tdHr  dHz\ 

=  i(ir-ir)  [5-23] 


dt 


dEz  c  ( Id 


dt 


e\r  dr''  ^  r  dip  J 


[5.2.5] 


dH, 


0 


dt 


=  c 


dEz  dEr 


dr 


_£M 

dz  J 


dHz 

dt 


1  dEr  Id 


=  c  — 


r  dp  r  dr 


[5.2.2] 


[5.2.4] 


[5.2.6] 


The  fields  can  generally  be  written  as 


f  EreikzZ~io)t  [5.2.7] 
E  =  Eoeikzz~io>t  J  Ecj)eikzZ~ia>t 
\E  e ikzz-iojt 


f  HreikzZ~io)t  [5.2.8] 
H  =  H0eikzZ~ib)t  -»  ]tf0ei/CzZ_iwt 

I  lj  pikzz—ia)t 


After  making  the  appropriate  substitutions  and  some  algebra  we  arrive  at  the  following  results 


Er 


i 

K2 


a)  dHz\ 
rc  dp  / 


ia)  dHz 
rK2c  dp 


E™ 


ikz  dEz 
k2  dr 


[5.2.9] 
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Hr 


i  /  eu>  dEz  dHz\ 
k2  V  rc  dtp  +  3r  / 


HTrE  = 


ikz  dHz 

k2  dr 
ieoo  dEz 

vk2c  dcp 


[5.2.10] 


f  _  ia)  dHz 
_  i  / kzdEz  codHz\  ^  J  K2C  gr 

^  K2\r  dcp  c  dr  J  |  „T m  _  dEz 

l  ^  rK2  dcp 


H  (j)  y 
v  K2 


i  t  eco  dE. 


z  kz  dHz > 


H, 


TE 


c  dr  r  dtp 


— » 


H, 


TM 


ikz  dHz 
rK2  dip 
ieoo  dEz 
k2c  dr 


[5.2.12] 


Using  these  equations  and  our  expression  for  Ez  ->  Ez  and  Hz  -*  Hz  we  can  now  write  the  fields 

in  the  r  —  direction  and  cp  —  direction,  let’s  note  that  all /m(JCjr)  ->  ]m  and  ( Khr )  -* 
for  simplicity 


iKz  V-1  10)  v-1 

Er(r  <R)  =  —  >  amJm  cos  (imp)  +  — j-  >  mbmJm  sin  (nup) 

k,  L-i  rKf  c  Z_i 

m  1  m 

,  ikl  v  .  , .  ,  in)  x-1  ,  , 

E\(r  >/?)  =  —  >  (a^H  +  +  amH  m)  cos (m0)  +  — j-  >  (/?£#«  +  pmHm)mcos(mcp) 

Ku  Z_j  rKuC/—i 


[5.2.13] 


[5.2.14] 


For  the  magnetic  field 


Hlr(r  <R)  = 


rtcfc 


^  '  Q-mJn 


msm(mcp )  + 


i/ci  v-1 

7^) 

Kin  ■  ^ 


sin  (rrup) 


[5.2.15] 
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,  ieh(jL iv-1  ,  ,  iki  v 

Hlr(r>R)  = - 2  /  «wm  +  amHm)  m  sin (m0)  +  —  >  (/?+//  +  +  /?mtf  m)  sin (m0) 

TKicZ-J  Kh  /—! 


[5.2.16] 


The  electric  fields  in  the  <p  —  direction 


[5.2.17] 


7  l«z  tU)  V-1 

EAr  </?)  =  -—  >  a^m  sin(?n0) - >  V/ m  sin(m0) 

T  KfrL-i  K;C  l—i 

1  m  m 

Etyir  >  R)  =  -ify  (amHm  +  a^H^msmimcp)  -  +  Pm^'m)  sin(m0)  [5.2.18] 

v  Kir  £—i  Khc  L-i 


The  magnetic  fields  in  the  <p  —  direction 


[5.2.19] 


,  iCjCn  v-1  ikz  v-1 

fi<*(r  <  /?)  = -  >  cos(m0)  H  2  /  hm/m^cos(m0) 

v  KjC  Z_i  rKT  /— I 

m  1  m 

fiqb(r  >  /?)  =  Y (a+//'+  +  a-fi'-)cos(?n0)  +  -^|V (p+H++p-H-)mcos(m(p)  [5.2.20] 

*  KhC  Z—i  rKt  Z-j 


5.3  Applying  the  boundary  conditions  at  r  =  R 

We  begin  by  using  the  continuity  of  the  tangential  components  of  the  electric  and  magnetic  fields 
at  r  =  to  obtain  an  expression  for  {am,  hm},  Z?^}  and  {cc^,  P^}.  In  the  z  —  direction 

Q-rrJm  ~  T  CCrnHrn  [5.3.1] 

KiJm  =  Pm^m  +  Pm^m 


[5.3.2] 
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and  in  the  0  —  direction 


n> 


<n 


KfR 


Q-mJrrjn  d"  ^m/  m  2  n  d-  d"  m  +  (3mH  m) 


<R 


KhC 


[5.3.3] 


€;(0  kl  ehOL)  kl 

—  arrJ'm  +  =  —  «H'+  +  «" //'"  )  +  (/M  +  /3mHm ) 

fCj  C  iiKj  « \.j^C  a 


[5.3.4] 


After  solving  Eqs.  (5.3.1)  and(5.3.2)  for  (am,  hm]  and  substituting  these  into  Eqs.(5.3.3)  and 
(5.3.4),  we  see  that  only  one  of  the  three  sets  of  coefficients  {am,  hm}  ,{<2^,  /?,},},  and  {<2,0  /?0}  is 
independent.  Explicitly,  the  linear  relationship 


[5.3.5] 


can  be  derived  from 


klm  (  1 


£i4w+ 

icWm  m 


LJr  + 


(Jin  1 


/c^m  ^  1 

”r”W 


U) 

c 


Jm  u- 

~rHr 

KiJm 


klm  (  1 

— U 
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7m  Ki 


h: 


[5.3.6] 
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where  each  of  the  four  sub-matrices  is  a  diagonal  matrix  with  its  elements  corresponding  to  the 
Bessel  function  combinations  evaluated  at  r  —  R.  In  this  cylindrically-symmetric  case,  the  S- 
matrix  can  be  formally  divided  into  four  (diagonal)  sub-matrices 


S  = 


Sn 

-S21 


S12 

$22. 


[5.3.7] 


The  components  5X1  and  S22  represent  polarization-preserving  TE-,  TM-reflection,  while  the 
components  $12’  $21  represent  polarization-mixing  coupling  of  TM  to  TE  waves.  Note  that  in 
the  cylindrical  geometry,  polarization-preserving  reflection  is  only  possible  when  either  m  =  0 
or  kz  —  0,  which  yields  detS12  =  detS21  =  0. 

Using  Eqns.  (5.3.1)  and  (5.3.2)  allows  one  to  calculate  the  amplitudes  {a,  b }  based  on  the 
amplitudes  of  {a~,  /?“}.  The  field  equations  will  provide  complete  information  about  the  field 
distribution  inside  the  unit  cell  once  the  parameters  klz,  and  { a~ ,  /?“}  are  known.  In  order  to 
obtain  the  additional  expressions  needed  to  solve  for  these  unknowns  the  periodicity  of  the  unit 
cell  must  now  be  investigated. 


5.4  Dispersion  of  the  Longitudinal  Mode 

Let’s  now  focus  on  the  problem  of  calculating  the  dispersion  of  the  mode.  This  reduces  to  the 
problem  of  calculating  a  relationship  between  internal  structure  of  the  unit  cell  and  the  set  of 
parameters  klz,  while  are  known.  For  the  square  unit  cell  geometry,  considered  in  this 

work,  the  latter  combination  will  only  contain  cylindrical  modes  with  m  =  0,4,8,...  The  field  of 
the  eigenmode  propagating  in  the  periodic  array  of  wires  should  satisfy  the  Bloch-periodicity 


condition 
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Here  we  enforce  periodicity  of  y  —components  of  electric  and  magnetic  fields.  Although 
this  condition  should  ideally  be  satisfied  for  all  values  of  the  y  coordinate  from  the  interval 

y  6  —  in  practice  it  suffices  to  enforce  the  Bloch-periodicity  condition  for  a  number  of 
fixed  points  {Xj,yf\  equal  to  the  number  of  m  values  in  Eqs.(5.4.2)  and  (5.4.3).  In  the 
calculations  we  assume  y,-  =  ——j,  with  Nm  being  the  number  of  m  tenns.  The  analysis 

2  Nm 

suggests  that  the  choice  of  the  exact  location  of  the  points  does  not  significantly  alter  the 
dispersion  of  the  mode,  derived  with  the  technique  described  below  (see  Fig  5.1). 
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Noting  from  symmetry  that  sin(0)  =  sin(7r  —  <p ),  cos (0)  =  —  cos(7r  —  0 ),  sin(m0)  = 
—  sin(m(7r  —  0)),  and  co s(m<p)  —  cos(m(n  —  0)),  it  can  be  shown  that  components  of  electric 
and  magnetic  field  possess  the  following  identities: 


Ey(x,  y)  =  Ey(—x,  y) 


Hy(x,y )  =  -Hy(-x,  y) 


Therefore,  Eq.(5.4.1)  becomes  equivalent  to 


[5.4.4] 

[5.4.5] 


f  0  0\  /a+\  (  0  6  \  (am\  _  f0\  [5.4.6] 

where  the  elements  of  the  sub-matrices  H±  are  evaluated  based  on  the  longitudinal  field 
equations  according  to  the  following  rules:  the  sub-matrices  H  and  H  represent  the  TM-  and  TE- 
polarized  magnetic  field  respectively;  the  superscript  of  the  expression  [+]  corresponds  to  the 
superscript  of  the  Hankel  function;  and  the  jm~th  element  of  the  sub-matrix  represent  the  y- 
component  of  the  magnetic  field  due  to  m-th  Hankel  function,  evaluated  at  the  point  (x, ,  yj  ]  = 

jj,  y;|.  With  the  help  of  the  S  matrix,  Eq.(5.4.6)  can  be  further  simplified  as 


VH+Sn  +  H+S21  +  H~ 


0  \  f^m\  _  /0\ 

iT+sl- 2  +  fns^2  +  fr)  \Pm)  ~  voy 


[5.4.7] 


Finally,  the  amplitudes  of  the  field  of  longitudinal  TM-polarized  wave  are  represented  as 
with  values  of  the  coefficients  a  given  by  non-trivial  solutions  of  the  linear  relationship 
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(Ky)a~  =  (H+Sn  +  H+S21  +  H a~  =  0  [5A81 

It  is  now  easy  to  see  that  the  dispersion  of  this  wave  is  given  by 

det|!Hy|=0  [5.4.9] 

To  verify  the  validity  of  Eq.(5.4.9),  the  dispersion  of  the  longitudinal  wave  corresponding  to  (i) 
m  —  0,4  and  (ii)  m  —  0,4,8  is  calculated  and  compared  with  the  derived  dispersions  to  numerical 
solutions  of  Maxwell’s  equations.  To  summarize  the  results  of  these  calculations  in  Fig.  5.1.  It 
is  clearly  seen  that  even  the  rough  approximation  with  m  —  0,4  yields  relatively  good  agreement 
with  numerical  solutions  of  Maxwell’s  equations.  Including  one  extra  harmonic  makes 
agreement  almost  perfect.  For  the  rest  of  this  work  we  use  the  results  corresponding  to  m  — 
0,4,8. 


Figure  5.1:  Propagation  constant  of  the  TM-polarized  waves  in  nano  wire  composite  as  a  function  of  wire 
permittivity.  Lines  and  dashes  represent  Eqn.  (5.4.9)  for  m=0,4,8  and  m=0,4  respectively;  symbols 
represent  the  numerical  solutions  to  Maxwell’s  equations. 

Now  that  the  longitudinal  dispersion  is  realized  the  remaining  coefficients  can  be  solved  and  an 
analytical  result  of  the  wave  profiles  can  be  produced.  These  results  are  in  excellent  agreement 
with  the  numerical  calculations  previously  obtained.  The  fields  generated  using  the  longitudinal 
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field  equations  in  comparison  with  the  FEM  field  profiles  show  striking  similarities  (Fig.  5.2) 
and  (Fig.  5.3). 


(a)  Re[El  (b)  Re[HJ 


-50  0  50  -50  0  50 

x  (nm)  x  (nm) 


Figure  5.2:  Analytical  electric  and  magnetic  fields  within  the  unit  cell  produced  using  the  previously 
derived  field  expressions,  (a,b)  corresponds  to  et  =  —  4  —  O.lt  and  (c,d)  to  =  —  9  —  O.lt. 

Notice  that  Since  (Ez)  =£  0,  this  verifies  that  the  solution  represents  the  longitudinal  wave  with 
dispersion  ez(/rz)  =  0. 

Fig.  5.3  demonstrates  the  excellent  agreement  between  the  numerical  and  analytical 
solutions  corresponding  to  a  three-term  series  m  —  0,4, 8,  and  clearly  demonstrates  the 
longitudinal  character  of  this  mode.  This  wave  is  strongly  dispersive  in  the  regime  -»  eh, 
corresponding  to  the  surface  plasmon  oscillations  on  the  metal-dielectric  interface.  On  the  other 
hand,  when  —  »  eh  (realized  at  mid-IR  and  lower  frequencies  for  noble  metals),  the 

wavevector  of  the  longitudinal  mode  approaches  n^oi/c,  and  its  transverse  counterpart 
approaches  the  light  line(see  [24]  and  Fig.2). 
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x  (nm)  x  (nm)  x  (nm) 

Figure  5.3 :  (a,b)  Dispersion  in  nanowire  composite  as  a  function  of  wire  permittivity.  Dashed  and  solid 
lines  represent  transverse  and  longitudinal  waves  k™9 and  klz,  eqn.  (5.4.9)  respectively,  symbols  represent 
numerical  solutions  to  Maxwell’s  equations;  for  —  et  »  €i,klz  ->  n^w/c  (dotted  line  in  (b)).  (c,d,e) 

electric  field  in  the  unit  cell;  surface  plots  and  arrows  represent  Ez  and  Exy  components,  respectively. 
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6  Solutions  at  Oblique  Angles 

The  next  step  is  to  obtain  the  dispersion  and  field  profiles  for  angles  other  than  incident.  To 
achieve  this  we  use  the  previous  derived  results  for  the  longitudinal  wave,  which  describes  the 
microscopic  field,  and  those  obtained  from  effective  medium  theory  (EMT),  the  macroscopic 
field.  Coupling  these  modes  together  through  the  use  of  the  coupled  oscillator  model  a  complete 
description  of  the  nonlocal  effective  permittivity  can  be  obtained. 

ez(kz)  =t(k2z-klz  )  — 

where  kz  is  the  wavevector  of  the  mode  in  the  nolocal  effective  medium  approximation,  klz  is 
the  wavevector  of  the  mode  composite  in  the  microscopic  theory,  and  is  the  factor  which  will 
be  determined  below. 

Eqn.  (6.1)  can  be  used  along  with  the  above  considerations  to  detennine  the  complete 
nonlocal  dispersion  for  propagation  at  an  angle  to  the  optical  axis.  For  simplicity,  we  consider 
the  case  ky  =  0,  kx  =£  0.  Using  egns.  (2.2.14)  and  (2.2.15)  where  the  propagation  of  the  TM- 

polarized  wave  can  be  described  by  £z(kz)  {kz  —  e™y  —  —e™ykx.  By  substituting  Eq.(6.1) 
in  this  equation,  we  obtain  the  following  relation 


mg 

X'V  c2) 


em9  ,  ,2 
cx,y  Oj 


2  kX- 


[6.2] 


Similar  to  other  nonlocal  materials,  nanowire  materials  support  two  TM-polarized  waves 
propagating  with  different  indices 
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The  final  free  parameter  of  the  model,  multiplicative  factor  ,  can  be  determined  by  requiring 
that  in  the  limit  of  small  kx  the  properties  of  one  of  the  two  TM-modes  follow  elliptical  or 
hyperbolic  dispersion  and  has  kz(kx)  =  const  dependence  that  is  observed  in  the  wire  media 
when  e ™9  >  0,  e™9  <  0,  and  e™9  «  — 1  respectively  [6,24,30,39,42,43].  The  relationship 


?  =  V 


ti  +  <Ti 


th  ~nt 


[6.4] 


adequately  describes  optics  of  wire  media  in  these  three  limits.  The  excellent  agreement 
between  the  predictions  of  Eq.  (6.3)  and  the  full- wave  numerical  solutions  of  Maxwell’s 
equations  are  shown  in  Fig.  6. 1 . 
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Figure  6.1:  (a,b)  Isofrequency  contours  of  TM-polarized  modes  in  nanowire  composites.  Solid  lines  and  symbols 
correspond  to  Eq.(6.3)  and  numerical  solutions  of  Maxwell  equations  respectively;  dashed  lines  represent  local  EMT. 

Cl) 

(c,d)  The  effecttive  medium  permittivity  of  nonlocal  nanowire  composite  for  kx  =  0  (c)  and  for  kx  ^  0  (d);  ex  , 
e ^  represent  two  solutions  of  Eq.(6.1);  cross  marks  ENZ  condition  (e™B  =  0). 


As  expected  the  isofrequency  of  the  “main”  TM-polarized  wave  resembles  and  ellipse  or 
hyperbola  which  for  small  values  of  kx  is  well  described  by  ema .  At  the  same  time,  the 
dependence  kz(kx)  for  the  “additional”  TM-wave  is  opposite  to  that  of  its  “main”  counterpart. 
The  z-component  of  the  permittivity  can  be  described  by  Eq.  (2.1.11)  when  kx  —  0,  but  for 
oblique  angles  exhibits  strong  spatial  dispersion. 

To  further  validate  our  analytical  approach  we  analyzed  the  optical  properties  of  the 
system  with  permittivity  of  the  host  material  eh  —  2.25.  Such  a  comparison  is  shown  in  Fig.  5.4. 
It  is  clearly  seen  that  the  developed  technique  is  in  excellent  agreement  with  numerical  solutions 
of  Maxwell  equations.  In  contrast,  local  EMT  adequately  predicts  response  of  the  nanowire 
system  only  outside  ENZ  region,  and  only  for  almost-normal  propagation  kx  «  oj/c. 
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Figure  6.2:  (a,b)  Isofrequency  contours  of  TM-polarized  modes  in  nanowire  composites  for  the  composite  with 
eh  =  2.25;  (c,d)  Isofrequency  contours  for  additional  TM-polarized  mode  for  the  same  system 


6.1  Wave  profiles  at  oblique  angles 

Now  that  the  origin  and  dispersion  of  the  modes  propagating  in  nanowire  systems  is  understood, 
let’s  focus  on  the  analysis  of  the  optical  properties  of  finite-size  wire  arrays.  Since  in  the  EMT 
approximation  the  fields  of  TE  and  TM-polarized  modes  are  orthogonal  to  each  other,  and  since 
propagation  of  TE  polarized  light  through  the  wire-based  system  is  only  affected  by  x,y- 
components  of  the  permittivity,  this  propagation  can  be  successfully  described  by  Eqs.(2.1.11) 
and  (2.1.12). 

Here,  we  focus  on  the  analysis  of  propagation  of  TM-polarized  light.  This  analysis  must  describe 
the  structure  of  electromagnetic  waves  propagating  in  the  system,  and  determine  what  additional 


46 


57 


boundary  conditions  needed  to  calculate  the  amplitudes  of  the  two  TM-polarized  modes  inside 
the  wire  system. 

Consistency  with  the  effective  medium  description  requires  that  the  unit-cell  averaged 
fields  satisfy  both  constituent  relations  £j  —  (e£)  )/(£)■)  and  relations  between  the  field 
components  of  the  plane  wave  ( eEz )  =  —kx/kz(eEx).  With  these  constraints,  we  start  with 
Eqs.(2.1.11)  and  (2.1.12)  and  the  field  solutions  derived  in  Ch.5.  One  can  determine  the 
parameters  ,  and  e™5  by  normalizing  (Elz)  —  (E™9)  =  (E™9)  =  1,  and  constructing  the 

fields  of  the  two  waves  propagating  in  the  wire  media  as  E(x,  y')eia>t~lkxX-lkzZ  with 

Ex{x,y)  =  E™9  [6.1.1] 

Ez(.x,y)  =  ym3  E™9  +ylElz  |z=0  [6- 1-2] 


Using  the  above  expressions  we  arrive  at 


ymg  _ 


fm9k 
cx,y  ^ x 


€ZkZ 


[6.1.3] 


em9  -  el 

cz  c 


yl  = 


pm9k  f 
c x,y  fcz 


£z  Ez 


[6.1.4] 


In  Eqs.(6.1.3)  and  (6.1.4)  ez  =  e(/cz)  which  is  given  by  eqn.  (6.1)  kz(kx)  by  eqn.  (6.3),  and 
el  =  ( e(x,y)Elz(x,y))/(Elz(x,y )>. 
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Eq.(6.1.2)  represents  a  transition  between  full-wave  solutions  of  Maxwell’s  equations  in 
the  nanowire  array  where  the  fields  oscillate  on  the  scale  of  the  individual  wires,  and  effective- 
medium  solutions  where  plane  waves  propagate  in  the  homogenized  material.  Since  our  model 
for  Em 9  assumes  quasistatic  limit,  Eq.(6.1.2)  is  technically  valid  in  the  limit  a  «  2n/kx,  A0. 
However,  Fig.  6.1  indicates  that  the  developed  fonnalism  also  provides  adequate  approximation 
for  optics  of  wire  systems  for  higher  values  of  kx. 
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7  Transmission  and  Reflection 

We  now  consider  the  problem  of  reflection/refraction  of  light  at  the  interface  of  two  (nonlocal) 
(meta-)  materials,  extending  the  well-developed  transfer-matrix  formalism  [56,57]  to  nonlocal 
composites. 

7.1  Transfer  Matrix  formalism 

The  typical  geometry  of  light  propagation  through  a  finite-thickness  slab  of  nanowire  material  is 
shown  in  Fig.7.1.  For  this  geometry,  Maxwell's  equations  require  continuity  of  (microscopic)  Ex 
and  Dz  fields.  Conventional  effective-medium  boundary  conditions  are  obtained  by  averaging 
these  relationships  across  the  unit  cell  [6,28,30,39,40].  However,  nanowire  media  support  two 
TM-polarized  modes  resulting  in  four  TM-polarized  waves  propagating  in  the  finite-thickness 
slab  (see  Fig.7.1).  Therefore,  calculation  of  the  amplitudes  of  these  waves  requires  additional 
boundary  conditions  (ABCs).  Different  forms  of  ABCs,  often  based  on  heuristic  arguments  have 
been  suggested  in  previous  works  [16,22].  Here  we  present  a  first-principles  approach  to  solve 
this  long-standing  problem. 


Figure  7.1:  Schematic  of  TMM  for  a  nanowire  composite.  Each  mode  is  labeled  to  clarify  notation  as  well  as 
interface  numbers. 
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Maxwell’s  equations  require  continuity  of  (microscopic)  Ex  and  Dz  and  the  effective- 
medium  boundary  conditions  and  is  obtained  by  averaging  these  relationships  across  the  unit 
cell.  Multiple,  linearly  independent  boundary  conditions  can  be  obtained  by  requiring  the 

continuity  of  En  —  (e2mna-Ex)  and  Dn  =  (eZninaDz)  with  different  integer  n. 

Continuity  of  E°  and  D°  yields  conventional  EMT;  continuity  of  D1  [and,  in  case  of 
contact  of  two  nonlocal  media,  of  E1]  represents  the  set  of  ABCs.  For  example,  assuming  that 
the  top  interface  is  located  at  z  =  z0,  the  boundary  conditions  for  this  interface  are  implemented 
as: 


chui,ie 


ikz^zo 

+  cxlExle~ 

ikz  1,1 

Zo 

-x 

C2,l 

E  °2,ieikz2’ 

iz° 

+ 

rj 

N  l 

E  he 

~lkz2 ,1 

z0 

1=1,2 

ifeZiiZ0 

—  c^D^e" 

Zo 

-I 

n°2leikz2 

lz  0 

~C2 ,1 

®he' 

~lkz2,\ 

,z0 

1=1,2 

ifczi.iZo 

—  cxl  Oi;1e" 

Zo 

-I 

ch 

M\,ieikz2 

lz0 

~  C2,l 

©2  ,le 

~lkz2,i 

|Z0 

1=1,2 


[7.1.1] 


In  the  expressions  above  the  double-subscript  represents  the  layer  number  in  the  system  and  the 
mode  number  within  the  layer,  while  the  “+”  superscript  represent  the  direction  of  the  wave 
propagation  (see  Fig.  7.1).  The  amplitudes  of  the  wave  propagating  in  the  system  represent  the 
amplitude  of  E°;  therefore,  the  amplitudes  of  Ex  are  symmetric  with  respect  to  change  of 
propagation  direction  kz  -»  —kz,  while  the  amplitudes  of  Ez,  Dz  are  anti-symmetric 

Eq.(7.1.1)  can  be  straightforwardly  used  as  a  building  block  to  implement  the  nonlocal 
analog  of  transfer-matrix  formalism,  calculating  transmission,  reflection,  and  field  distributions 
across  an  arbitrary  planar  stack  containing  nanowire  composites. 

As  before,  two  sets  of  numerical  calculations  were  used,  commercial  finite-element 
solver  and  in-house  RCWA  code.  Since  both  sets  of  calculations  utilized  full  3D  geometry  of  the 
unit  cell  of  finite-length  wire  array,  and  since  solutions  of  Maxwell’s  equations  in  3D  geometry 
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are  known  to  be  challenging,  we  analyzed  the  stability  of  results  of  our  calculations.  To  assess 
the  accuracy  of  RCWA  calculations,  the  number  of  Fourier  hannonics  in  representation  of  spatial 
profile  of  permittivity  was  varied.  Similarly,  the  mesh  size  in  finite-element  calculations  is  varied 
and  the  stability  of  the  solutions  with  respect  to  minor  modification  of  the  geometry  are  checked. 
In  particular,  the  cap  parts  of  the  wires  were  changed  from  the  planar  to  ellipsoidal  shape  [see 
Fig. 7. 2].  The  results  of  these  studies  are  summarized  in  Fig.7.2.  Note  that  according  to  finite- 
element  calculations,  transmission  and  reflection  profiles  strongly  depend  on  the  geometry  for 


6t  —  —3.  Analysis  of  the  field  profiles  in  this  spectral  range  suggests  that  these  minima  are 
artifacts  of  meshing  routines  implemented  in  the  FEM  solver.  To  mitigate  these  artifacts,  we 
show  the  average  of  these  results. 


h=2nm  h=4nm 


RCWAh=0nm  <>• 
FEM  h=0nm  <■ 
FEM  h=2nm  »  ▼ 

FEM  h=2nm  a  a 


-10  -9  -8  -7  -6 


Figure  7.2:  (a)  Schematic  of  different  nanowire  geometry  used  in  finite-element  numerical  solutions  of  Maxwell 
equations;  inset  shows  to-scale  geometry;  (b-d)  Transmission  and  reflection  of  light  through  a  parallel  slab  of 
nanowire  media,  suspended  in  air  with  nanowires  with  /m(e;)  =  —0.1  at  various  angles  for  different  nanowire 
geometry  and  solution  method. 
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Transmission  and  reflection  of  the  nanowire  metamaterial  are  now  compared  for  full- 
vectorial  numerical  solutions  of  Maxwell's  equations  predictions  of  conventional  (local)  effective 
medium  theory,  and  predictions  of  the  nonlocal  EMT  developed  in  this  work  in  Fig. 7. 3.  It  is 
seen  the  smaller  the  loss  and  the  larger  the  angle  of  incidence  the  more  important  it  becomes  to 
take  into  account  the  nonlocal  optics  of  nanowire  composites.  Interestingly,  nonlocal  response 
strongly  affects  optical  response  of  the  wire  metamaterials  across  the  broad  range  of  the  effective 
permittivities.  This  effect  most  clearly  seen  in  reflection,  but  is  also  visible  in  transmission, 
especially  in  ENZ  regime. 


Figure  7.3:  Transmission  and  reflection  of  light  through  a  parallel  slab  of  nanowire  media,  suspended  in  air  with 
Im{ei')  =  —0.1  (a,b)  and  /m(e;)  =  —0.25  (c,d).  (a,c):  local  TMM  calculations,  (b,d):  nonlocal  EMT  developed  here 
(lines)  and  numerical  solutions  of  Maxwell  equations  (symbols);  Solid  lines  and  filled  symbols  represent  reflection, 
dashed  lines  and  empty  symbols  -  transmission. 


To  further  validate  the  developed  approach,  we  calculated  transmission  and  reflection  for 
a  planar  slab  of  nano  wire  material  with  permittivity  of  the  host  eh  —  2.25.  A  comparison 


52 


63 


between  the  numerical  solutions  of  Maxwell’s  equations,  predictions  of  local  EMT,  and  nonlocal 
EMT  is  shown  in  Fig. 7. 4.  We  again  note  that  the  nonlocal  EMT  provides  a  substantially  more 
accurate  description  of  optics  of  nanowire  systems  than  its  local  counterpart. 


e.  e. 

I  l 


Figure  7.4:  Transmission  and  reflection  of  light  through  a  parallel  slab  of  nano  wire  material  suspended  in  a  substrate 


with  eh  =  2.25  and  Im{e{)  —  —0.1.  (a)  Local  calculations  and  (b)  nonlocal  EMT  (lines)  and  numerical  solutions  to 


Maxwell’s  equations  (symbols). 


To  assess  the  robustness  of  the  developed  formalism,  we  compared  transmission  and  reflection 
of  light  by  nanowire  arrays  calculated  with  the  developed  nonlocal  EMT  and  with  numerical 
solutions  of  Maxwell’s  equations. 
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8  Conclusion 

To  conclude,  an  approach  to  describe  the  optics  governing  nonlocal  wire  metamaterials  across 
the  whole  optical  spectrum  was  presented.  This  fonnalism  demonstrates  excellent  agreement 
with  the  results  of  numerical  solutions  to  Maxwell's  equations  and  is  essential  in  the 
development  of  an  adequate  description  of  optics  in  wire  arrays.  This  facilitates  understanding 
the  true  density  of  photonic  states  through  the  limits  of  resolution  in  these  systems. 

The  formalism  can  be  expanded  to  calculate  the  field  inside  the  structures  with  more 
complicated  unit  cells,  that  include  multi-shelled  wires,  coax-cable-like  systems,  etc.  [52,53]  and 
can  be  further  extended  to  calculate  fields  inside  the  systems  with  a  non-circular  cross-section  of 
wires,  in  which  case  the  S-matix  ceases  to  be  block-diagonal  [54,55].  A  mechanism  is  also 
provided  to  explain  the  deviation  for  the  dispersion  of  the  waves  in  nanowire  materials  from  the 
prediction  of  local  effective  medium  technique.  In  particular,  this  deviation  places  fundamental 
limits  on  subwavelengh  light  manipulation,  and  on  the  increase  of  the  local  density  of  photonics 


states  [21]. 
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9  Future  Work 

The  developed  analytical  model  that  adequately  describes  the  steady-state  single-frequency 
response  for  nanowire  systems  can  be  straightforwardly  be  adapted  to  wire  structures  with  non¬ 
square  unit  cells  as  well  as  to  wire-like  composites,  e.g  multi-shelled  wires,  and  coax-cable-like 
systems.  The  fonnalism  can  also  be  used  as  a  platform  for  analysis  of  complex  light  interaction 
with  anisotropic  plasmonic  composites. 

Our  next  major  task  is  to  form  an  understanding  of  the  light  emission  in  nanowire 
metamaterials.  It  is  well  known  that  light  emission  by  dyes,  quantum  dots,  and  other  quantum 
mechanical  systems  can  be  substantially  altered  by  the  surrounding  materials.  Due  to  the 
correlation  between  light  emission  and  the  density  of  states,  any  modulation  of  the  local  density 
of  states  is  known  to  lead  to  enhancement  (or  quenching)  of  fluorescence  rates  of  emitters  in 
proximity  to  metal  nanostructures  [58-61]. 

Emission  properties  are  further  modified  in  hyperbolic  systems  that  in  effective  medium 
approximation  possess  singular  density  of  states  over  a  significant  frequency  range  [62,63]. 
However,  while  it  is  clear  that  nonlocality  will  significantly  affect  emission  in  metamaterials 
with  large  (or  small)  values  of  effective  parameters,  the  very  definition  of  density  of  optical 
states  in  nonlocal  systems  is  not  completely  understood.  Future  work  is  to  develop  an 
understanding  of  interaction  between  nonlocal  metamaterials  and  light-emitting  systems. 
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